Longitudinal and transverse meson correlators 
in the deconfined phase from the lattice 

Gert Aarts*, Chris Allton 1 *, Justin Foley t , Simon Hands* and Seyong 

Kim** 

* Department of Physics, Swansea University, Swansea SA2 8PP, United Kingdom 
^ Department of Physics, Carnegie Mellon University, Pittsburgh, PA 15213, USA 
£N| | ** Department of Physics, Sejong University, Seoul 143-747, Korea 

> 

o 

/■"■n ■ Abstract. It has long been known that QCD undergoes a deconfining phase transition at high 

temperature. One of the consequent features of this new, quark-gluon phase is that hadrons become 
, unbounded. In this talk meson correlation functions at non-zero momentum are studied in the 

deconfined phase using the Maximum Entropy Method. 
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QCD is well-known to be a strongly interacting and confining theory under "normal" 
conditions. However, it has long been theorised that at large energy scales (in temper- 
ature, T, or baryonic chemical potential, /l) it undergoes a transition to a deconfined 
(quark-gluon plasma) phase. This has been experimentally observed at CERN [1] and at 
the RHIC experiment at Brookhaven [2]. The exact nature of this transition is undergo- 
ing intense theoretical and experimental investigation. In fact, even the general features 
of the QCD phase diagram are still being mapped out (see Fig.l). As an example of this 
uncertainty, the Particle Data Book [3] does not contain a single reference to the de- 
confined phase of QCD ! Properties of QCD at T = fl = are measured experimentally 
and calculated theoretically often to accuracy at the percent level or below. The same is 
certainly not true of the deconfined phase, where properties typically have much larger 
(~ 20%) errors associated with them, if they are known at all. 

Physics at T = /i = is characterised by quantities such as hadronic masses and 
transition matrix elements, whereas the relevant quantities in the deconfined phase 
correspond to those of plasma physics: pressure, entropy, susceptibilities, and response 
functions. Of particular interest in this work are transport coefficients. 

From general arguments based on the asymptotically free nature of QCD, the naive 
expectation is that quarks and gluons would be virtually free in the deconfined phase. 
However, at RHIC, a (relatively) strongly interacting phase was found with "almost 
instantaneous" equilibration and a small ratio of viscosity, r\, to entropy density, s (both 
signs of "strong" coupling). These are characteristics of a so-called "Perfect Fluid". 
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FIGURE 1. Illustrative phase diagram of QCD. 



There is some theoretical basis for such small values of rj/s in the plasma phase. A 
calculation relying on the correspondence between jV = 4 supersymmetric Yang-Mills 
theory and superstring theory in A dSs x S 5 space-time predicts a lower bound for r]/s 
[4,5], 

7]/s>— for N c , g 2 N c — >• °°. (1) 

While this theory is not QCD, it gives a scale for this ratio in a strongly interacting 
theory [5]. 

Transport coefficients are essential properties of the plasma phase and can be derived 
theoretically from current-current spectral functions. We list some of them here. 

Shear viscosity, r\ , which is obtained from off-diagonal energy-momentum correlators 

in the zero energy limit [6]. 
Bulk viscosity, £ , which is obtained from diagonal energy-momentum correlators in 

the zero energy limit [7]. 
Electrical conductivity, o, and Diffusivity, D, which are obtained from the energy 

dependence of vector meson spectral functions [8]. 

The aim of this talk is to make progress in the understanding of the deconfining 
mechanism by analysing mesonic systems simulated by a lattice calculation at T ^ (but 
with fi = 0). Naively, it is expected that the two (valence) quarks in these mesons should 
become unbound at the deconfining temperature, T c . In practice this appears true for light 
quark states only, and charmonium states appear to remain bound until temperatures 
higher than T c [9, 10]. Specifically we aim to determine the diffusivity through a lattice 
simulation of correlation functions of vector currents at non-zero momentum. This work 
extends our previous studies at zero momentum [8] where we calculated the electrical 
conductivity. A full version of this work is in preparation [11]. 



LATTICE BACKGROUND 



The conventional approach to studying hadronic quantities with the lattice technique is 
via the imaginary-time dependence of Euclidean correlation functions of operators with 
well-defined quantum numbers. In the confined phase, each bound state (in the tower of 
states with those quantum numbers) contributes ~ e~ Eit to the correlation function, G{t), 
where E\ is the state's energy. At large times, t, the lowest state dominates, so fitting G{t) 
to an exponential form can trivially (in theory) determine Eq. 

In the deconfined phase, the situation is more subtle. Unbound states no longer 
contribute pure exponential terms to G(t), and it is more appropriate to introduce the 
spectral function, p(co,p), 

f 00 dco 
G(t,p) = j p(co,p) K(t, co) — , (2) 

where the (lattice) kernel is defined 

coshK>-i/(2r))l 

sinh[ft)/(2r)J 

and we have allowed for a momentum dependence in the correlation function G and 
therefore in p. As usual, the temperature, T is the inverse temporal length, 1 / (aN t ). 

In the confined case, each non-decaying state i contributes a delta function, <5(ft) — 
Ei), to p(co). A decaying state would have a spectral feature of finite width, and an 
unbound state would correspond to a continuous spectrum. By studying the temperature 
dependence of spectral functions the transition from the bound to deconfined phases can 
be observed. 

As well as containing information on the stability or otherwise of hadronic states, 
p(ft),/>) also can be used to extract transport coefficients (as described in Sec.l) and 
hydrodynamic structure. 

However, despite the importance of the spectral function, and its simple definition in 
terms of the correlation function in eq.(2), it is notoriously difficult to extract. This is be- 
cause it is an example of an ill-posed problem: there are (in general) more ft) data points 
in p(ft>) than there are t data points in the correlation function G{t). The method which 
has now become fairly standard to overcome this problem is the Maximum Entropy 
Method (MEM) which is based on Bayesian statistics (for a review, see [12]). MEM is 
a very standard technique in fields which require image reconstruction/deconvolution, 
such as astronomy, crystallography, and the analysis of atomic/molecular spectra. 

This work is an extension to our earlier work [8, 13, 14] to non-zero momentum. 
Specifically, we apply MEM to calculate p(ft),p) for mesonic correlation functions on 
lattice data with parameters summarised in Table 1 . We used the quenched approxima- 
tion with the standard Wilson gluonic and staggered fermionic actions, full details are 
given in [8]. Twisted boundary conditions [15, 16] were used to allow a finer resolution 
in momentum space for G(t,p) (and therefore p(ft),p)). In all 21 different momenta 
combinations were studied (ranging up to \p\ ~ 10/L) of which 17 are non-degenerate. 

We follow our earlier work [8] where a singularity in K((D,t) as ft) — > was corrected 
by a simple redefinition of K — > G)K/{2T). This allows us to obtain reliable p(ft),j?) 
estimates, even in the ft) — > limit. 



TABLE 1. Lattice parameters used in the simulation. 







Cold 


Hot 


Spatial Volume 


Nf x N, 


48 3 x 24 


64 3 x 24 


Lattice spacings 


a- 1 


-4GeV 


- 10 GeV 


T 


l/{aN t ) 


r~ 160MeV-0.62r c 


T ~420MeV~ 1.5r e 


Statistics 


Nc f8 


- 100 


- 100 




FIGURE 2. An example of the staggered correlation function showing the MEM analysis. 



RESULTS 

Our fermionic action uses staggered quarks so therefore the hadronic correlators are 
a mixture of states with opposite parity, see Fig. 2. These can be decomposed into the 
following spectral representation, 

GM)= /o ^^ w )[p( £0 ^)-(- 1 )'p( £0 ^)]- W 

This means that to recover the physical spectral function, the even and odd timeslices 
must be treated separately and then combined, 

P = ^(p eV6n + p 0dd ). (5) 

Since there are only N t = 24 time slices in our lattices, the MEM analysis of the even 
and odd timeslices includes six timeslices (allowing for time reversal symmetry). This 
motivates the use of anisotropic lattices in future studies [9]. 
The electrical conductivity, defined 
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FIGURE 3. Longitudinal and transverse vector correlation functions normalised by the average corre- 
lation function, G^ ve = Gj^ + 2Gfs for various momenta, p, as a function of time. Data points for each 
momenta (for a given time) are offset horizontally for clarity; from left to right they are \p\L = 0, 2, n, 2n, 

was found using this method [8], where pa is the spectral function for the spatial 
component of the vector correlator. The work presented here extends [8] by including 
non-zero momenta. In Fig. 3, we plot the longitudinal and transverse vector correlation 
function versus time for various momenta. As can be seen, there is a distinct difference 
between the longitudinal and transverse correlators and a clear systematic effect as the 
momenta increase. 

The diffusivity, D, can be obtained, in principle, from the momentum dependency of 
the longitudinal vector spectral function in the light quark mass limit (see e.g. [17]). 
In Fig. 4, the spectral functions for the longitudinal vector case are shown for both the 
ma = 0.01 and 0.05 quark masses. As can be seen, there is a clear non-zero intercept in 
the case of the 0.01 mass which is absent in the 0.05 case. It is this non-zero intercept in 
the 0.01 case which led us to determine the conductivity in [8]. Our future plans are to 
study this momentum dependency with the aim of independently determining D [11]. 
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